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Abstract 

I present results from an approach that extends the Eliashberg theory by systematic expansion in the vertex function; 
an essential extension at large phonon frequencies, even for weak coupling. In order to deal with computationally 
expensive double sums over momenta, a dynamical cluster approximation (DCA) approach is used to incorporate 
momentum dependence into the Eliashberg equations. First, I consider the effects of introducing partial momentum 
dependence on the standard Eliashberg theory using a quasi-local approximation; which I use to demonstrate that it 
is essential to include corrections beyond the standard theory when investigating d-wave states. Using the extended 
theory with vertex corrections, I compute electron and phonon spectral functions. A kink in the electronic dispersion is 
found in the normal state along the major symmetry directions, similar to that found in photo-emission from cuprates. 
The phonon spectral function shows that for weak coupling WX < coo, the dispersion for phonons has weak momentum 
dependence, with consequences for the theory of optical phonon mediated d-wave superconductivity, which is shown 
to be 2nd order in A. In particular, examination of the order parameter vs. filling shows that vertex corrections lead 
to d-wave superconductivity mediated via simple optical phonons. I map out the order parameters in detail, showing 
that there is significant induced anisotropy in the superconducting pairing in quasi-2D systems. [PUBLISHED AS: 
JOURNAL OF PHYSICS AND CHEMISTRY OF SOLIDS, VOL. 69, 2982-2985 (2008)] 
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1. Introduction 

Angle-resolved photo-emission spectroscopy 
(ARPES) directly probes the dispersion of electrons; 
identifying a kink associated with the active optical 
phonon in the cuprates [1]. Neutron scattering has 
shown an anomalous change in the phonon spec- 
trum at the superconducting transition, indicating 
an interesting role for phonons in the cuprates [2], 
Estimates of the magnitude of the electron-phonon 
coupling, and the energy of the phonon mode put 
the problem outside the limited region of applica- 
bility for BCS theory, so schemes to cope with a 
wider range of parameters need to be developed. 
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Moreover, since there a node in the superconduct- 
ing gap consistent with d-wave symmetry [3], any 
theory implicating phonons as the mechanism must 
be able to explain the unconventional order. 

Electron-phonon (c-ph) interactions can be de- 
scribed by the following generic model, 



kq 



]>> q (& q 6 q + 1/2) 



(1) 



Here, h = 1 and ek = —2t(cos(ak x ) + cos(afcj,)) — 
2t±cos(ck z ) with t — 0.25 and t± = 0.01 (rep- 
resenting a quasi-2D system and taming the van- 
Hove singularities). The Holstein model is a spe- 
cial case where <? q = go and a-> q = luq, for which 
A = Sft/O^o^O is the dimensionless e-ph coupling 
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(W is the half-band width). The momentum- 
independent phonon dispersion approximates an 
optical phonon, and the momentum-independent 
e-ph coupling corresponds to a completely local 
interaction. A Coulomb pseudo-potential, = 
-W^ c J2kqa( c ticT c -qa) c k<jC-]<<7 + h.c., may be added 
(i.e. Hubbard model at the Hartree-Fock level). 
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Fig. 1. (a) Overview of the diagrammatic representations 
of the self-energy used in this work. £ is the electron, and 
II the phonon self energy. The row labeled ME contains 
the vertex neglected and VC the corrected diagrams, (b) 
Schematic of the DCA (quasi-local) formalism for cluster 
size Nq = 4. Within the squares S(k) i— > S(K) where K 
are represented by solid circles. Equivalently, A(k) i— » A(K). 
A partial density of states (DOS) is associated with each 
sub-zone. DCA enables computations in the thermodynamic 
limit, via an equivalent cluster impurity problem. 

Some treatments of the ARPES kink used a stan- 
dard Eliashberg theory with anisotropic coupling 
and an assumed d-wave order parameter for fitting 
purposes [4,5]. The extended Eliashberg approach 
that I have developed differs significantly from pre- 
vious work in the sense that it is fully self-consistent, 
with no prior assumptions about the form of the or- 
der parameter (i.e. any ci-wave order arises directly 
from the self-consistency of the theory) [6,7]. 

2. Eliashberg theory and momentum: Failure 
for optical-phonon mediated eZ-wave states 

In this section, I discuss how quasi-local Eliash- 
berg equations can be constructed, and demonstrate 
how the lowest-order Fock diagram does not con- 
tribute to the gap equations. The calculation is in- 
tended to be demonstrative, and shows that stan- 
dard Eliashberg equations fail for an optical phonon 
mediated cZ-wave state. 

The Eliashberg equations are computed from the 
Fock contribution to the self-energy: 

E(k, iw n ) = -T^ J j0^glf 3 G(k- q,w n ') 

xf 3 K(q, iuj n - icj n i) (2) 



and equate S with the approximate form, 
S w (1 - Z k )iw„f + A k (iw„)fi + xt~3 



(3) 



Normally the Eliashberg equations are computed 
using a local (momentum-independent) approxima- 
tion. However, to examine ci-wave states, the lo- 
cal approximation is not sufficient. Substituting the 
Green function G(k, zw„) _1 = iuJ n To — Ck^3 — S in 
eq. 2 (£k = ft — ^ are Pauli matrices and [i is the 
chemical potential), and treating the gap function 
as constant in sub-zones of k-space in the manner of 
DCA, one obtains a quasi-local approximation for 
the Eliashberg equations, 
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When gap functions have d-wave symmetry, A^o = 
— Ao i7r and Ao.o = A^.^ = 0, and the partial DOS 
2?(0, 7r) = V(tt, 0) has the symmetry of the lattice, it 
is immediately clear that the 0(A) contribution to 
the quasi- local Eliashberg equations vanishes, and 
2nd order terms must be considered. When there 
is some weak momentum dependence in the e-ph 
coupling and phonon frequency, this cancellation 
will not be exact, but the contribution from 2nd- 
order terms can still be expected to be the largest 
in the perturbation series. Thus a fully extended set 
of Eliashberg equations including vertex corrections 
needs to be constructed. Currently, it is not possi- 
ble to construct a similar approximation to eq. 4 for 
the vertex corrected theory. However a numerical 
approach using the full DCA equations is available 
to enlighten some aspects of the e-ph problem. 

The rest of this article concerns results from an 
extended Eliashberg theory with momentum depen- 
dence and vertex corrections determined numeri- 
cally using DCA. No form for the order parameter 
is assumed in advance and full forms for the Green 
functions and self-energies are used. Full details of 
the translation of the diagrams in Fig. 1(a) may be 
found in ref . [6] . 

3. Normal state spectroscopy 

Figure 2 shows image plots of the electron spec- 
tral functions [8]. Vertex corrections are most pro- 
nounced for ujo > WX, as expected from analysis 
of the parameter space [9]. In particular, the vertex 
corrected theory shows a kink at the energy scale 
relating to the phonon frequency measured at the 
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(a) A(o>), ME, a> =0.4W, X=0.3 (b) A(o>), VC, co =0.4W, X=0.3 




Fig. 2. (Colour online) Electron spectral functions comparing VC and ME approximations when cluster size Nc = 64. The 
effects of vertex corrections are most pronounced for wq > WX, as expected from analysis of the parameter space [9]. Horizontal 
dashed lines indicate the renormalised phonon frequency at the (it, 0) zone-edge as computed within the self-consistent theory 
(i.e. the one which would be measured in a neutron experiment). In particular, only the vertex corrected theory shows a kink 
at the energy scale relating to the phonon frequency measured at the (?r, 0) point [8] as reported in Ref. [1]. 



(-7T, 0) point (i.e. renormalised zone edge phonon fre- 
quency) in agreement with the ARPES measure- 
ments in Ref. [1] . 

Phonon spectral functions are shown in Fig. 3. 
The phonon renormalisation is extremely small for 
c^o = 0.4W. In contrast, modes at the (ir, ir) point 
in the Brillouin zone are strongly softened for luq — 
0.2W when vertex corrections are not taken into ac- 
count. The second order terms act against the soft- 
ening and reinforce the theory. A small reduction in 
spectral weight is seen mid way along the T-X and 
T-W lines when luq > WX. This might be a precur- 
sor to the strong softening seen in Ref. [2] on pass- 
ing through the transition temperature, but further 
analysis of the superconducting state is necessary to 
be sure. 

4. Evolution of the order parameter 

d-wave superconductivity is one of the remarkable 
results from the theory extended with vertex cor- 
rections. Figure 4 shows the anomalous density as 
the chemical potential is varied first upwards from 
half-filling and then back from large filling showing 
that two solutions are stable. A small fj, c = 1.2 is 
applied. A = 0.6 and loq/W — 0.4. Starting with 
a half-filled solution and increasing, then decreas- 
ing the chemical potential, first d-wave and then (on 
decreasing from high electron density) s-wave solu- 
tions to the self-consistent equations are found (as 
detailed in the caption of Fig. 4) . Highly anisotropic 
solutions form as the Fermi-surface approaches the 
van-Hove points. In the d-wave state, the first order 
terms are not the leading order of the perturbation 
theory in A. It was argued in Ref. [7] that this was in 
part due to the nearly flat phonon spectrum, so that 
the phonon propagator D(k, u>) w D{w) and the 



1st order contribution to the anomalous self-energy 
WX J2q D(q)F(k — q) is small (F is the anomalous 
Green function) . The anomalous contribution to the 
pseudopotential term is (f>c = Wfio ^(k)- Since 
F is modulated in the d-wave state, then the con- 
tribution of the pseudopotential term is zero. In the 
s-wave state, F is not modulated, and the result- 
ing finite contribution destroys the s-wave order for 
large enough nc, thus selecting a dominant d-wave 
contribution. The phonon spectra in Fig. 3 demon- 
strate that so long as luq > WX, then this condi- 
tion is met, and 2nd order terms are essential. 2nd 
order terms remain the leading order correction to 
the anomalous self-energy, even for the larger A used 
here. It is clear that the 2nd order terms are similar 
for both repulsive and attractive models, and thus 
it is perhaps not surprising that both attractive and 
repulsive models show a d-wave state. 

5. Summary and Outlook 

I have described spectroscopic results and order 
parameters from an extended Eliashberg-style the- 
ory with vertex and momentum corrections. The 
extended theory is essential for the investigation 
of d-wave superconductivity mediated by optical 
phonons, since the vertex corrections are the low- 
est order terms in the perturbation expansion in 
electron-phonon coupling A. Indeed, d-wave super- 
conductivity is immediately predicted with such an 
extension. 

There are two remaining areas which require fur- 
ther investigation. The first involves a better treat- 
ment of the Coulomb repulsion, which is currently 
only treated as a simple Hartree-Fock pseudopoten- 
tial. In particular, the ability to treat large repul- 
sion and antiferromagnetism is essential for high- 
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o((o), ME, X=0.3, to W = 0.4 




Fig. 3. (Colour online) Phonon spectral functions. The phonon renormalisation is moderate for ujq = OAW, Nc = 64. White 
points show the values computed from uj 2 — o>q(1 — Re[n(u;)]) = for clarification, (a) Modes at the W point in the Brillouin 
zone are almost completely softened in the theory without vertex corrections (a split mode is indicated in panel (b)). Vertex 
corrections act against that softening and reinforce the theory (c,d). There is an indication of reduction of spectral weight half 
way between the high symmetry points when ljq > WX, perhaps with similar origin to that in the neutron results [2] 




Fig. 4. Anomalous density as /i is varied first upwards from half-filling and then back from large filling showing two types of 
solutions to the self-consistent equations, fic = 1.2, A = 0.6, Nq = 4 and loq/W = 0.4. (a) Around half-filling (n = 1), there is 
only one solution with d-wave order. This persists to n = 1.175 where the order parameter spontaneously changes to s-wave 
order (b,c). The reduction in d-wave magnitude is reminiscent of a quantum critical transition. When decreasing fi from dense 
filling, the order parameter is a conventional s-wave (d), which persists below n = 1.175, i.e. there are 2 solutions (e). Below 
n = 1.05 there is only one solution. Close to n = 1.05, the s-wave solution has an extended-s character, as can be seen in panel 
(0- 



Tc materials. Treatment of momentum dependent 
e-ph coupling would also be of interest. However, 
the overriding challenge at this stage is to convert 
the numerical results into a coherent analytic theory 
with vertex corrections. 
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